Quantum theory of tunneling magnetoresistance in 
GaMn As / GaAs / GaMn As heterostructures 

Alireza SaffarzadetQ and Ali A. Shokri 

Department of Physics, Tehran Payame Noor University, 
Fallahpour St., Nejatollahi St., Tehran, Iran 
(Dated: February 6, 2008) 

Abstract 

Using a quantum theory including spin-splitting effect in diluted magnetic semiconductors, we 
study the dependence of tunneling magnetoresistance (TMR) on barrier thickness, temperature 
and applied voltage in GaMnAs/GaAs/GaMnAs heterostructures. TMR ratios more than 65% 
are obtained at zero temperature, when one GaAs monolayer (~ 0.565 nm) is used as a tunnel 
barrier. It is also shown that the TMR ratio decreases rapidly with increasing the barrier thickness 
and applied voltage, however at high voltages and low thicknesses, the TMR first increases and 
then decreases. Our model calculations well explain the main features of the recent experimental 
observations. 
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I. INTRODUCTION 



In the past few years, tunneling magnetoresistance (TMR) in magnetic tunnel junctions 
(MTJs) has attracted much attention due to the possibility of its application in magnetic 
memories, magnetic field sensors and quantum computing devices In a common MTJ, 

the TMR between two ferromagnetic electrodes separated by an insulator layer depends on 
the spin polarization of the ferromagnetic electrodes and spin dependent tunneling rates 
For such structures, much efforts have been made to increase the TMR ratio, 
which can be defined by the relative change in the current densities. Recently, a large TMR 
ratio as much as 60% has been reported in Co-Fe/Al-O/Co-Fe junctions [6] which opened a 
new era of technical applications of the tunnel junctions. 



Using ferromagnetic semiconductors (FMSs 



vestigated in single 0, 0, 0, 13] and double 
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such as EuS, the TMR has also been in- 
magnetic barrier junctions. In 



these structures, the FMSs which act as spin filters are used as tunnel barriers. Therefore, 
the MTJs using these materials, are able to produce highly spin-polarized current and high 
values for TMR. However, the Cure temperature in most of FMSs is much lower than room 
temperature jl5||. 

On the other hand, coexistence of ferromagnetism and semiconducting properties in di- 
luted magnetic semiconductors (DMSs) has opened the prospect of developing devices which 
make it possible to combine the information processing and data storage in the same mate- 
rial jlfil . \v\ . In the magnetic semiconductors, the exchange interaction between the itinerant 
carriers in the semiconducting band and the localized spins on magnetic impurity ions leads 
to spectacular magneto-optical effects, such as giant Faraday rotation or Zeeman splitting. 
Among different types of DMSs, Gai_ x Mn x As is one of the most suitable materials for use 
in spintronic devices. In this kind of materials, a fraction of the group-Ill sublattice is 
replaced at random by Mn ions which act as acceptors. According to the strong Hund cou- 
pling, the five 3d electrons in a Mn 2+ ion give rise to a localized magnetic moment. The 
DMSs based on III-V semiconductors, like Gai-aMn^ As compounds, exhibit ferromagnetism 
with Curie temperature as high as 110 K for Mn doping of around x=0.05. The origin of 
ferromagnetism in the GaAs-based DMSs can be demonstrated through the p-d exchange 
coupling between the itinerant holes in the valence band of GaAs and the spin of mag- 
netic impurities [3, 12, Q- In spite of low Curie temperature of Gai_ x Mn x As, the use of 
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Gai-^Mn^As-based III-V heterostructures as MTJs presents several advantages. The possi- 
ble advantages include not only simple integration with existing semiconductor technology, 
but also relatively straightforward fabrication of high-quality epitaxial structures, easily con- 
trolled physical and electronic structure, and the integration of quantum heterostructures 
which is easier than in any other materials system |ig| | . 

The spontaneous magnetization and the feasibility of preparing GaAs-based DMS 
are stimulated much attention for studying the spin-polarized transport phenomena in 
Gai_3.Mn3.As/ (GaAs or AlAs) heterostructures. Recently, Chiba et al. 2l( observed a 
TMR ratio of 5.5% at 20 K in a GaMnAs tunnel junction. Also, Tanaka and Higo [3] 
reported TMR ratios more than 70% at 8 K in GaMnAs/ AlAs/GaMnAs heterostructures. 
Moreover, recent analyses of temperature dependence of current-voltage characteristics of 
GaMnAs/GaAs/p-GaAs p-i-p diodes have shown that GaAs intermediary layer which is a 
nonmagnetic semiconductor (NMS), acts as a barrier (87-140 meV) for holes injected from 
GaMnAs 2_l. in this regard, a large TMR ratio of 290% has been observed at 0.39 K 
in GaMnAs/ GaAs/GaMnAs heterostructures, around zero applied bias (2^). Such a high 
TMR ratio is a result of high spin polarization of DMS electrodes and the high quality of 
sample and interfaces between GaMnAs and GaAs layer. Interlayer exchange interaction is 
another effect which has also been studied in Gao.96Mno.04As/GaAs superlattices, where an 
oscillatory behavior with varying thickness of GaAs spacer is predicted j^j]. 

From a theoretical point of view, the TMR has also been investigated in GaM- 
nAs/ AlAs/GaMnAs heterostructures. Tanaka and Higo [22] calculated the dependence of 
TMR on AlAs thickness using the tight-binding theory including that the parallel wave 
vector of carriers is conserved in tunneling process. The results showed that the TMR ra- 

n 

tio decreases rapidly with increasing barrier thickness. Also, Tao et al. 26] studied the 
AlAs thickness dependence of TMR using the transfer matrix approach. Their treatment, 
however is somewhat questionable, because they did not apply the boundary conditions for 
derivatives of wave functions correctly. When the materials are different, it is the normal 
mass flux that must be continuous. 

In the present work, we study theoretically the spin-polarized transport in GaM- 
nAs/ GaAs/ GaMnAs tunnel junctions by considering the effect of spontaneous magnetization 
of GaMnAs layers. The temperature dependence of the spin splitting energy is calculated 
for DMSs within the mean-field approximation. Using the transfer matrix method and the 
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nearly-free-carrier approximation we will study the TMR and spin polarization of tunneling 
carriers in terms of the barrier thickness and bias voltage at all temperatures. We assume 
that the carrier wave vector parallel to the interfaces and the carrier spin are conserved in 
the tunneling process through the whole system. 

The paper is organized as follows. In section 2, the model is described and then the 
spin polarization and the TMR for DMS/NMS/DMS structures are formulated. Next, the 
temperature dependence of spin splitting energy is obtained for DMS, within the mean- 
field framework. In section 3, numerical results obtained for a typical tunnel junction are 
discussed. The results of this work are summarized in section 4. 



II. DESCRIPTION OF MODEL AND FORMALISM 

In order to investigate the spin-dependent transport properties in MTJs based on DMS 
materials, we consider two semi-infinite DMS electrodes separated by a NMS layer, which 
acts as a barrier, in the presence of DC applied bias V a . For simplicity, we assume that the 
two DMS electrodes are made of the same materials and all of the interfaces are flat. In the 
framework of the parabolic-band effective mass approximation, the longitudinal part of the 
one-hole Hamiltonian can be written as 

where m* (j=l-3) is the hole effective mass in the jth layer, and 



Uj{x) 



0, x < , 

E F + (p- eV a x/d, < x < d , (2) 
— eV a , x > d , 



where d is the barrier thickness, Ep is the Fermi energy in the left electrode, measured from 
the middle point between the edges of the two spin subbands, and <fi is the barrier height 
measured from the Fermi level. — h^ F • a is the internal exchange energy where h^ F is the 
molecular field in the jth DMS electrode and a is the conventional Pauli spin operator. 

The Schrodinger equation for a biased barrier layer can be simplified by a coordinate 
transformation whose solution is the linear combination of the Airy function Ai[p(x)] and 
its complement Bi[p(x)] [221 ■ Considering all three regions of the DMS/NMS/DMS junction, 
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the eigenfunctions of the Hamiltonian (0) with eigenvalue E x have the following forms: 



A Xa e iklaX + B la e~ iklaX , x < , 
A 2a Ai[p(x)] + B 2a Bi[p(x)], < x < d, 
A 3a e ik3 ° x + B 3a e- ik3 * x , x > d, 



(3) 



where, 



k la = J2m\(E x + h a)/h 



k 



(T) 



(4) 
(5) 



v 3CT - y^m^ + eV a + Th a)/h , 

are the hole wave vectors along the x axis. Here, a are the hole spin components ±1 (or 
t, I), ho = |h^ F | and T = +1(— 1) for parallel (antiparallel) alignment of the magnetizations. 
The coefficients Aj a and Bj a are constants to be determined from the boundary conditions, 
while 



p{x) 



d 



x 



eV a X 



(E F + (f>-E x - -eV a 



d 



with 



A 



-h 2 d 



2m* 2 eV a 



1/3 



(6) 



(7) 



Although the transverse momentum kii is omitted from the above notations, the summa- 
tion over k|| is carried out in our calculations. 

Upon applying the boundary conditions such that the wave functions and their first 
derivatives are matched at each interface point Xj , i.e. ipj^ixj) = i>j+i,<r( x j) an d 
(m*)~ l [di/jj !(T (xj)/dx} = (m* +l )~ 1 {dil)j + i a (xj)/dx], we obtain a matrix formula that con- 
nects the coefficients A\ a and B\ a with the coefficients A 3a and B 3a |8[. Since there is no 
reflection in region 3, the coefficient -830- in Eq. (jH)) is zero and the transmission coefficient 
of the spin a hole, which is defined as the ratio of the transmitted flux to the incident flux 
can be written as 



TP(E x ,V a 



(r) 



Am\m* 3 k la k 3a 



ak la kfj + 8m{m* 3 ) 2 + ((3m* 3 k la - ^m\k K 3a 



(r)x2 



where the following abbreviations are used 

a = Ai[p(0)]Bi[p(d)] - Bi[p(0)]Ai[p(d)], 

1 
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Am 2 
1 



Am 



-{Ai[p(0)]Bi'[p(d)]-Bi[p(0)]Ai'[p(d)]}, 
{Ai>(0)]Bi[p(d)]-Bi'[p(0)]Ai[p(d)]}, 



(\ m * 2 y 



■{Ai'[p(0)]BiV(d)]-BiV(0)]Ai>(d)]}- 



(8) 

(9) 
(10) 

(11) 
(12) 



Here, Ai'[p(x)] and Bi'[p(x)] are the first derivatives of the Airy functions. 



A. Current densities and TMR 



The spin-dependent current density for a MTJ with a given applied bias V a and at 
temperature T can be calculated within the nearly-free- hole approximation [28]: 



» = f 1F(E„ K)ln { 1 + 7/%- r . iE. , (13) 

47r 2 ^ Jej [ 1 + exp l(E F - E x - eV a )/k B T\ J 

and, at T=0 K, 

jP(Va) 



eV a TP(E x ,V a )dE x + (E F -E x )TP(E x ,V a )dE x , 

J El JE F -eVa 

(14) 

where ks is the Boltzmann constant and Eq is the lowest possible energy that will allow 
transmission and is given by Eq = max{— ho, —eV a — Tho} for spin-up holes and Eq = 
ho for spin-down ones. It is clear that the tunnel current is modulated by the magnetic 
configurations of the both DMS electrodes. 

The degree of spin polarization for the tunnel current is defined by P = (J^ — Ji)/(J^ + Ji), 
where (J;) is the spin-up (spin-down) current density. For the present structure, we obtain 
this quantity when the magnetizations of two DMS electrodes are in parallel alignment. For 
studying the TMR, the spin currents in the both parallel and antiparallel alignments are 
calculated. In this case, the TMR can be defined as 



( J t p + J?) - ( J T ap + Jf p ) 

' I'M n = ^ — )a P + jl P — l — , ( i ■ ) ) 



where Jl ' , corresponds to the current density in the parallel (antiparallel) alignment of 
the magnetizations in the FM electrodes, for a spin-up (f ) or spin-down (j) hole. 



B. Temperature dependence of spin-splitting energy 

In the following, we develop the formalism for calculating of temperature dependence of 
the spin splitting energy in the DMS electrodes using the mean-field theory. At the first 
step, we consider the magnetic Hamiltonian for one itinerant hole with spin s located at 
r and one Mn ion with spin S located at R. Both of these magnetic subsystems can be 
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described as 

H h = Jpd^s-S/^r-Rj) , (16) 
i 

and 

# Mn = j P dE s *- s <K r *- R ) • ( 17 ) 

i 

Here, J p d is the hole-Mn exchange coupling constant, fi B is the Bohr magneton, g h and gjjn 
are the Lande g-factor for holes and Mn ions, respectively. In the mean-field approximation, 
H h — ► H& F = g h fi B s ■ h h and H Mn — > = gun^S ■ h Mn , where h h and h Mn are 

respectively the effective magnetic fields acting upon holes and magnetic impurities, and 
can be given as 



and 



h h = —J pd N Mn (S z ) , (18) 

gh^B 



Here, Nun = 4x/a[] is the density of Mn ions, with a 3 , being the unit cell volume. M h = 
J2i( s i) 3( r i — R-) is the magnetization density of the hole subsystem, which is assumed to be 
uniform within the length scale of the magnetic interactions, so the magnetic response of 
the the Fermi sea holes to the effective field h h , is given by 

M h = g ht i B s 2 D h (E F )h h , (20) 

where the density of states of the hole gas with effective mass m* and the hole concentration 
p is D h (Ep) = (37T 2 ) _2 / 3 (3m*/^ 2 )p 1 / 3 . On the other hand, the magnetic response of the 
impurity spin to the effective field h Mn , is given by 

M Mn = g M n»BN Mn (S z ) = gun^NunSBs ( ^^— \ , (21) 

where 

Baix) = coth(^±^) - (-L) coth(^) , (22) 

is the Brillouin function. Therefore, within the spirit of a mean-field framework the magne- 
tization of Mn subsystem is then given by 

M Mn = g M ^ B N Mn SB s (^M h ) , (23) 
which should be determined self-consistency with the hole magnetization 

M h = JpdL>h(EF) M Mn . (24) 

2gUn^B 
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Also, the ferromagnetic transition temperature Tq can be obtained using the expansion for 
the Brillouin function B,s(x), when x C 1. In this regard, Tq is found as 

2 S(S+1)D(E F ) 

ic = {JpdS) iVMn • (25) 

6Kb 

Now, the temperature dependence of spin splitting energy in the jth DMS electrode can be 
obtained as Aj = 2\hf F \, where, the molecular field is hf F = (J pd /2g Mn fi B )Mf n . 



III. NUMERICAL RESULTS AND DISCUSSIONS 

Numerical calculations have been carried out to investigate the effects of barrier 
thickness, temperature and applied voltage on spin-dependent transport in a typical 
Gai-^Mn^As/GaAs/Gai-^Mn^As tunnel structure. We have chosen Gai-^Mn^As and 
GaAs because of the same crystal structures and lattice constants |22|. The relevant pa- 



rameters for Gai_ x Mn x As electrodes are chosen as ao = 5.65 A, E-p=0.2 eV 



23, J pd =0-15 



eV-nm 3 |30j, 5=5/2, s=l/2, g h = # M n=2, p = 4.9 x KPcm" 3 31] and T c = 110 K for a 
sample with x = 0.05 The suitable parameter for the barrier height of GaAs is 0=0.1 
eV I2J. 

It is important to note that in the valance band of semiconductors, there are heavy hole 
and light hole bands which are degenerate at the top of the valance band. Therefore, it 
seems to be reasonable that we use of one value for the hole mass at the structures with low 
Mn concentration. Hence, the effective mass of all carriers are taken as m*=0.16 (m^ 
is the free-hole mass), which is consistent with the experimental parameters have been used 
in the present work. 

In Fig. 1 we have presented the TMR as a function of the applied voltage for several 
temperatures when d = 0.565 nm. As it is evident from the figure, at all temperatures 
the TMR decreases monotonically as V a increases. Also, it is seen that, the bias voltage 
values, where the TMR ratio reaches half its maximum, varies from 55 to 90 mV depending 
on temperature. This value is much smaller than that of metal-based single barrier MTJ 
structures (300-400 mV) flHQ- 

At very low voltages, for both parallel and antiparallel 
alignments, the tunnel current densities vary linearly, therefore the TMR behaves similar to 
a flat curve. However, with increasing the applied voltage, due to the modification of the 
barrier shape, a nearly parabolic dependence of currents on the voltage appears, and this 
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modifies the TMR curve. With increasing the temperature, the spontaneous spin splitting 
in the valence band of the DMS layers, decreases. In this case, the difference between 
current densities in both parallel and antiparallel configurations is also reduced. Thus, in 
such temperatures one can expect low values for TMR ratio, as it is clear in the figure. The 
results show that the discrepancy between the TMR(100 K) and TMR(0 K) decreases with 
increasing the bias voltage. Therefore, the applied bias has a strong influence on the TMR. 
The relative behavior of the curves in Fig. 1, is qualitatively in good agreement with the 
recent experimental results j24(. 

The barrier thickness dependence of TMR at different applied biases has been shown in 
Fig. 2. For V a =5 and 80 mV, with increasing the barrier thickness, the TMR decreases 
rapidly from its maximum value and always remains positive for various thicknesses of the 
GaAs layer. However, for V a =120 and 150 mV, it is found that the TMR becomes negative 
in a certain range of the barrier thickness, as it is shown in the inset. The origin of this effect 
is related to the enhancement of the applied bias and can be observed, when the applied 
bias is more than the barrier height (100 meV). In fact, at high- voltage range (eV a > <p) a 
quantum well will appear at GaAs/GaMnAs interface; thus, the carriers at the Fermi level 
of the left electrode tunnel through the GaAs barrier and the quantum well, into the right 
electrode. In this regard, the carriers can cause standing waves in the barrier, leading to a 
reduction and a sign reversal of the TMR, as it is shown in the inset of Fig. 2. Also, with 
further increasing of the applied voltage, the TMR exhibits a peak at low thicknesses, which 
has been observed previously in GaMnAs/AlAs/GaMnAs heterostructures 

The temperature dependence of the TMR and spin polarization is also calculated in the 
present MTJ. The results are shown in Fig. 3 at V a =5 mV when d = 0.565 nm. From 
the figure, it is clear that, both the TMR and spin polarization decrease with increasing 
temperature, and vanish at 110 K, which is the Curie temperature of the Gao.95Mno.05As 
layers. One can see that, the highest value of TMR (about 70%) and spin polarization (about 
37%) can be obtained at zero temperature. The inset shows the normalized magnetization 
of Gao.95Mno.05As layers as a function of temperature. The reduction of TMR and spin 
polarization arises as a result of the dependence of spin splitting energy on temperature. The 



comparison of the present results and those obtained with experiment |2J] shows somewhat 
different behavior at 10-30 K. Such discrepancy between our theory and the experiment 
can be probably due to the effects of magnetic anisotropy at the GaMnAs/GaAs interfaces 



9 



and also spin-flip scattering at low temperature in the GaMnAs electrodes. The mean-field 
theory does not take into account the effects mentioned here. 

Our final comment addresses the role of Rashba spin-orbit interaction in the present 



tunnel structure [3^, |35[ . The Rashba effect produces a spin-flip scattering in the interfaces, 
which affect on the transmission of each spin channel and can not be solved separately. 
Therefore, a spin mixture occurs and the currents become less spin polarized. However, the 
Rashba effect is very small. In fact, with assuming the Rashba interaction in the magnetic 
tunnel junctions, it has been shown that, the spin- dependent transmission coefficients change 
very small in the both Parallel and antiparallel configurations Hence, the variation of 
spin polarization and TMR due to this effect is nearly small. Therefore, one can omit this 
effect in such tunnel heterostructure. That is, the Rashba effect does not have a dominant 
influence on the spin polarization and TMR. 



IV. SUMMARY AND CONCLUSIONS 



On the basis of the effective mass approximation model including the tem- 
perature dependence of spin splitting, we studied the spin-polarized tunneling in 
Gao.95Mno.05As/GaAs/Gao.95Mno.05As heterostructure. The numerical results showed that, 
the TMR strongly depends on GaAs barrier thickness, applied voltage and temperature. At 
high voltages, with increasing the barrier thickness, a peak and sign reversal observed in 
the TMR curves. Also, we found that, both the TMR and spin polarization decrease with 
increasing the temperature. The results obtained in this work show the main features of 
the experiments and can be a base of development for spin dependent tunneling electronic 
devices. 
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FIG. 1: Dependence of the TMR on the applied voltage at different temperatures. 



13 




[%] a i/Mi 

FIG. 2: Variations of the TMR as a function of the barrier thickness at different temperatures. 
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FIG. 3: Dependence of the TMR (left) and spin polarization (right) as a function of tempera- 
ture. The inset shows temperature dependence of the normalized magnetization of Gao.95Mno.05As 
electrodes. 
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